Draft version November 24, 2011 

Preprint typeset using I^TeJC style emulateapj v. 11/10/09 



PUMPING THE ECCENTRICITY OF EXOPLANETS BY TIDAL EFFECT 

Alexandre CM. Correia 

Department of Physics, I3N, University of Aveiro, Campus de Santiago, 3810-193 Aveiro, Portugal; and 
Astronomie et Systemes Dynamiques, IMCCE-CNRS UMR8028, 77 Av. Denfert-Rochereau, 75014 Paris, France 



GWENAEL BOUE 

Centro de Astrofisica, Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal; and 
Astronomie et Systemes Dynamiques, IMCCE-CNRS UMR8028, 77 Av. Denfert-Rochereau, 75014 Paris, France 

AND 

Jacques Laskar 

Astronomie et Systemes Dynamiques, IMCCE-CNRS UMR8028, 77 Av. Denfert-Rochereau, 75014 Paris, France 

Draft version November 24, 2011 

ABSTRACT 

Planets close to their host stars are believed to undergo significant tidal interactions, leading to a 
progressive damping of the orbital eccentricity. Here we show that, when the orbit of the planet is 
excited by an outer companion, tidal effects combined with gravitational interactions may give rise 
to a secular increasing drift on the eccentricity. As long as this secular drift counterbalances the 
damping effect, the eccentricity can increase to high values. This mechanism may explain why some 
of the moderate close-in exoplanets are observed with substantial eccentricity values. 
Subject headings: celestial mechanics — planetary systems — planets and satellites: general 



1. INTRODUCTION 

Close-in exoplanets, as for Mercury, Venus and the ma- 
jority of the natural satellites in the Solar System, are 
supposed to undergo significant tidal interactions, re- 
sulting that their spins and orbits are slowly modified. 
The ultimate stage for tidal evolution is the synchroniza- 
tion of the spin and the circularization of the orbit (e.g. 
Correia 2009). The spin evolves in a shorter time-scale 
than the orbit, so long-term studies on the tidal evolu- 
tion of exoplanets usually assume that their rotation is 
synchronously locked, and therefore limit the evolution 
to the orbits. However, these two kinds of evolution can- 
not be dissociated because the total angular momentum 
must be conserved. Synchronous rotation can only occur 
when the eccentricity is very close to zero. Otherwise, the 
rotation rate tends to be locked with the orbital speed 
at the periapsis, because tidal effects are stronger when 
the two bodies are closer to each other. In addition, in 
presence of a companion body, the eccentricity under- 
goes oscillations (e.g. Mardling 2007), and the rotation 
rate of the planet shows variations that follow the eccen- 
tricity (Correia & Laskar 2004). As a consequence, some 
unexpected behaviors can be observed, such as a secu- 
lar increase of the eccentricity. In this Letter we provide 
a simple averaged model for the orbital and spin evolu- 
tion of an exoplanet with a companion (Sect. 2), and ap- 
ply it to the HD 117618 planetary system (Sect. 3). We 
then give an explanation for the eccentricity pumping 
(Sect. 4), and derive some conclusions (Sect. 5). 

2. THE MODEL 

We consider here a system consisting of a central star 
of mass mo, an inner planet of mass mi, and an outer 
companion of mass m2. We use Jacobi canonical coordi- 
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nates, with ri being the position of mi relative to mo, 
and r2 the position of m2 relative to the center of mass of 
mi and mo- We further assume that |ri| <C |r2|, that the 
system is coplanar, and that the obliquity of the planet 
is zero. The inner planet is considered an oblate ellipsoid 
with gravity field coefficients given by J2, rotating about 
the axis of maximal inertia, with rotation rate c<;, such 
that (e.g. Lambeck 1988) 



J2 = h 



3Gmi 



(1) 



G is the gravitational constant, R is the radius of the 
planet, and k2 is the second Love number for potential. 

Since we are interested in the secular behavior of 
the system, we average the motion equations over the 
mean anomalies of both orbits. The averaged poten- 
tial, quadrupole- level for the spin (e.g. Correia & Laskar 
2010a), octopole- level for the orbits (e.g. Lee & Peale 
2003; Laskar & Boue 2010), and with general relativity 
corrections (e.g. Touma et al. 2009) is given by: 



[7 = -Co(l-e2)-i/2-Ci(l-e2)- 



■3/2 



-Co 



+ Q 



(l-ei)3/2 ^ (l-ei)5/2 



C0SU7 , (2) 



where 
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ai is the semi-major axis (that can also be expressed 
using the mean motion n^), is the eccentricity, and 
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w = voi— W2 is the difference between the longitudes of 
the periastron, Wi. We also have /3i = momi/(mo + mi), 
and 1^2 (mo + mi)m2/(mo + mi + m2). 

The contributions to the orbits are easily obtained us- 
ing the Lagrange planetary equations (e.g. Murray & 
Dermott 1999): 



e. = ^l^, * = -^?^. (5) 



Thus, 



62(1 + 3/46^)^1^ . 
-^31 _ ^2)5/2 ^1^^ ' 

ei(l + 3/4e^) . 

62 = ^^32 ^^_^2)2 ^ ' (^) 



and 
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where x = uj/ni^ and the constant frequencies 
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Notice that the variations in ei and 62 (Eqs. 6, 7) are re- 
lated by the conservation of the total angular momentum 
(after dividing by min2a2): 

n2 \a2j mi n2 \a2j ^ mi 

(15) 

where ^ is a structure coefficient. The conservative sys- 
tem (Eq. 2) can thus be reduced to one degree of freedom. 

In our model, we additionally consider tidal dissipation 
raised by the central star on the inner planet. The dissi- 
pation of the mechanical energy of tides in the planet's 



interior is responsible for a time delay At between the 
initial perturbation and the maximal deformation. As 
the rheology of planets is badly known, the exact depen- 
dence of At on the tidal frequency is unknown. Several 
models exist (for a review see Correia et al. 2003; Efroim- 
sky & Williams 2009), but for simplicity we adopt here 
a model with constant At, which can be made linear 
(Singer 1968; Mignard 1979). The contributions to the 
equations of motion are given by (e.g. Correia 2009): 



ni 

Oi 

ai 



-K(/i(ei)x-/2(ei)) 



2K' (/2(ei)x-/3(ei)) , 



ei=9KM^/4(ei)^-/5(ei))ei, 



where 



K = ni 



3k2 mo/3i 
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K' 



K mi 
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(17) 



(18) 



(19) 



Q-i = niAt, and /i(e) = (1 + ie^ + 3eV8)/(l - e^f'^, 
/2(e) = (l + 15eV2 + 45eV8 + 5eVl6)/(l-e2)6, /3(e) = 
(1 + 31eV2 + 255e^/8 + 185eVl6 + 25eV64)/(l - 6^)15/2, 
/4(e) = (l + 3eV2 + eV8)/(l-e2)5, /5(e) = (l + 15eV4+ 
15eV8 + 5eV64)/(l-e2)i3/2. 

We neglect the effect of tides over the longitude of the 
periastron, as well as the flatenning of the central star. 
Their effect is only to add a small supplementary fre- 
quency to wi , similar to the contributions from the gen- 
eral relativity (for a complete model see Correia et al. 
2011). 

Under the effect of tides alone, the equilibrium rotation 
rate, obtained when w = 0, is attained for (Eq. 16): 



- = /(ei) = ^ = l + 6ef + 0(et). 

ni /i(ei) 



(20) 



Usually K' <C so tidal effects modify the rotation 
rate much faster than the orbit. It is thus tempting to 
replace the equilibrium rotation in expressions (17) and 
(18). With this simplification, one obtains always nega- 
tive contributions for di and ei (Correia 2009), 



Oi 

ai 



-7K' U{ei)el , 



ei = --K'h{ei){l-el)ei 



(21) 



(22) 



with /6(e) = (l+45eVl4+8e4+685eV224+255eV448+ 
25eiVl792)(l - e^)-^^^^/{l + Se^ + 3eV8). Thus, the 
semi-major axis and the eccentricity can only decrease 
until the orbit of the planet becomes circular (Fig. la). 
However, planet-planet interactions can produce eccen- 
tricity oscillations with a period shorter, or compara- 
ble to the damping timescale of the spin. In that case, 
the expression (20) is not satisfied and multi-planetary 
systems may show non-intuitive eccentricity evolutions, 
such as eccentricity pumping of the inner orbit (ei in- 
creases while 62 decreases). 
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TABLE 1 

Single planetary systems with 0.1 < ai < 0.3 and ei > 0.3. 



Star 


ai 


ei 


mi 


mo 


Age 


r 


(name) 


(AU) 








(Gyr) 


(Gyr) 


HD 108147 


0.102 


0.53 


0.26 


1.19 


2.0 


0.01 


CoRoT-10 


0.105 


0.53 


2.75 


0.89 


3.0 


0.24 


HD 33283 


0.145 


0.48 


0.33 


1.24 


3.2 


0.34 


HD 17156 


0.163 


0.68 


3.19 


1.28 


3.4 


0.44 


HIP 57050 


0.164 


0.31 


0.30 


0.34 




39.4 


HD 117618 


0.176 


0.42 


0.18 


1.05 


3.9 


2.06 


HD 45652 


0.228 


0.38 


0.47 


0.83 




93.3 


HD 90156 


0.250 


0.31 


0.06 


0.84 


4.4 


35.8 


HD 37605 


0.260 


0.74 


2.84 


0.80 


10.7 


10.6 


HD 3651 


0.284 


0.63 


0.20 


0.79 


5.1 


15.5 



Note. — Data is taken from http://exoplanet.eu/ 



2miaf (l-ef )^ 
At21A;2Gm2j?5 



with At = 10^ s (Correia & Laskar 2010b). 



3. APPLICATION TO EXOPLANETS 

As an illustration of the eccentricity pumping, we apply 
our model to the HD 117618 system. This Sun-like star 
(mo ~ Mq) has been reported to host a single Saturn- 
like planet on a eccentric orbit (Butler et al. 2006). The 
residuals of the best fitted solution to the observational 
data are 5.5 m/s, so we assume that any additional com- 
panion with a doppler shift semi-amplitude smaller than 
this value is presently undetected, that is, any planet 
with 7712 < 0.2 Mj and a2 > 1.4 AU. 

In our simulations we adopt for the observed planet the 
same geophysical parameters as for Saturn, R = Rsati 
k2 = 1/2, ^ = 1/5, and a dissipation time lag At = 200 s 
(which is equivalent to Q ~ 3 x 10^). Since the semi- 
major axis of the planet undergoes tidal dissipation, its 
value was certainly larger when the system formed. We 
then adopt ai = 0.25 AU as initial value for all sim- 
ulations. The initial eccentricity is chosen such that 
ei ~ 0.4 when ai = 0.175, the present observed val- 
ues (Tab. 1). We further assume initial w = 180°, and 
27r/cj = 50 day, that quickly evolves near the equilibrium 
rotation (Eq. 20). 

In absence of a companion, the eccentricity and the 
semi-major axis are damped following an exponential de- 
cay (Correia & Laskar 2010b), and the present configu- 
ration is attained after 1 Gyr (Fig. la). At present the 
observed eccentricity would be around 0.1, and we still 
needed to explain the high initial value near 0.7. 

We now add a companion to the system with m2 = 
0.2 Mj, a2 = 1.8 AU, and 62 = 0.4, and set ei = 0.3. 
At first, we only consider dissipation in the spin (Eq.l6) 
and neglect its effect on the orbit (Eqs. 17, 18), in order 
to highlight the eccentricity pumping (Fig. lb). We then 
clearly observe this effect, the eccentricity of the inner 
planet rising up to 0.7. We also observe that the ec- 
centricity of the outer planet is simultaneously damped, 
because of the conservation of the total angular momen- 
tum (Eq. 15). 

Orbital and spin evolution cannot be dissociated, so 
we then integrate the full system (Fig. Ic). We observe 
that the initial behavior of the system is identical to the 
situation without dissipation on the orbit (Fig. lb). How- 
ever, as the eccentricity increases, the inner planet comes 
closer to the star at periastron, and tidal effects on the 
orbit become stronger. As a consequence, the semi-major 
axis decreases and the damping effect on the eccentricity 



(Eq. 18) overrides the pumping drift. The system ulti- 
mately evolves into a circular orbit. The present config- 
uration is attained around 4 Gyr of evolution, which is 
compatible with the present estimated age of the star. 
The pumping effect is then responsible for a delay in the 
final evolution of planetary systems and may explain the 
high values observed for some of them (Tab. 1). 

Finally, we repeat the integration of the full system, 
but with a smaller-mass companion m2 = 0.15 Mj at 
a2 = 1.4 AU. The companion eccentricity is still 62 = 0.4, 
but the inner planet now begins with ei = 0.15 (Fig. Id). 
The initial evolution is still similar to the previous sim- 
ulation (Fig. Ic), except that the eccentricity oscillations 
of both planets are higher, because the orbits are closer. 
Since the companion mass is smaller, its eccentricity also 
decreases more than before, and reaches zero around 
3 Gyr. At this stage, the angle w stops circulating, and 
begins librating around 180°. The two orbits are then 
tightly coupled and evolve together, showing an identi- 
cal behavior to close-in planets with a few days of orbital 
period (e.g. Mardling 2007; Laskar et al. 2011). As a 
consequence, the evolution time-scale is much longer, al- 
lowing the inner planet to maintain high eccentricity for 
longer periods of time (Fig. Id). The present eccentricity 
is only observed after 6 Gyr. 



4. ECCENTRICITY PUMPING 

In order to understand the unexpected behavior of the 
eccentricity during the initial stages of the evolution, we 
can perform some simplifications in the equations of mo- 
tion without loss of generality (Sect. 2). We can neglect 
tidal effects on orbital quantities (Eqs. 17, 18), which is 
justified since K' <^ K (Eq.l9). The only contribution 
of tides is then on the rotation rate (Eq.l6). The semi- 
major axis and the mean motion are thus constant, and 
the eccentricity only varies due to the gravitational per- 
turbations (Eq. 6). In addition, we linearize the set of 
equations of motion in the vicinity of the averaged val- 
ues of X, ei, and 62. Let x = + fe, where xq is the 
solution of (20), ei = eio + ^ei, and 62 = 620 + ^^2. In the 
following, 5e2 is expressed as a function of 5x and 5ei us- 
ing the conservation of the angular momentum (Eq.l5). 
Then, the set of equations of motion (6, 8, 16) reduces 
to: 

5ei = -Asinw , (23) 



Sx = —iyx^x-\- Ve^^l 5 



with 



i^31 



e2o(l + 3/4e2o)7l-e2 
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(l-eio)^/^ 



(24) 
(25) 

(26) 



(1 



°io) 



vixl 



(1 



2 Y 



-10 



-10 



(l-eio)3/2 



^22 
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(1 
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Fig. 1. — Long-term evolution of the HD 117618 system in different situations. We show the semi-major axis ai in AU (blue), and 
the eccentricities ei (red) and 62 (green), (a) without a companion; (b) with a 1712 = 0.2 Mj companion at 02 = 1.8 AU, but without 
dissipation on the orbit; (c) same as (b), but with a full model; (d) with a m2 = 0.15 Mj companion at 02 = 1.4 AU, with a full model. 
The dashed line gives the slope of the linear approximation (Eq. 37). 
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where /{(e) = 15(e + 36^2 + eV8)/(l - e^^^^, and 
/^(e) = 3(9e + 656^2 + 125 / 8 + 5e78)/(l - . We 
neglected the octupole terms since <C i^2i, the contri- 
butions from (5e2, and assumed that e^o / 0. 

At first order, the precession of the periastron is con- 
stant w and the eccentricity is simply given from 
expression (23) as 



5ei = Aecos{gt + vuq) , 



(32) 



where Ae = A/g. That is, the eccentricity ei presents 
periodic variations around an equilibrium value eio, with 
amplitude Ae and frequency g. Since gxSx^geSei <C g, 
the above solution for the eccentricity can be adopted 



as the zeroth order solution of the system of equations 
(23—25). With this approximation, the equation of mo- 
tion of 6x (25) becomes that of a driven harmonic oscil- 
lator whose the steady state solution is 



Sx = Ax cos{gt + wq 



(33) 



with Ax = UeAe/y^ -\- g'^^ and sin0 = gj ^Jv^ + g^ - 
The rotation rate thus presents an oscillation identical 
to the eccentricity (Eq.32), but with smaller amplitude 
and delayed by an angle ^ (see Correia 2011). Using the 
above expression in equation (24) and integrating, gives 
for the periastron: 

w = gt-\-wo-\-^Axsin{gt-\-wo — (j))-\-—Aes'm{gt-\-wo) . 

^ ^ (34) 

Finally, substituting in expression (23) and using the ap- 
proximation g^Ax^geAe <C g gives 



5ei ~ — A sm.{gt + 

9xA 
9 



— — Ae sm{gt + vuq) cos{gt + wq) 
Ax s'm{gt -\- wq — (/)) cos{gt + ujq) , (35) 
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or, combining the two products of periodic functions, 
Sei = —A sm{gt + wq) 1 — Ax sin(2^t + 2wo — 0) 

-^Aesm{2gt -\- 2zuo) + ^Axsin(/) . (36) 

The two middle terms in the above equation can be ne- 
glected since they are periodic and have a very small 
amplitude {gxAx^ geAe <C g). However, the last term in 
sin(/) is constant and it adds a small drift to the eccen- 
tricity. 



< ei >= 



t . 



(37) 



The drift is maximized for g ^ it vanishes for weak 
dissipation (Ax — ^ 0), but also for strong dissipation 
(0 0). Note that the phase lag between the ec- 
centricity (Eq. 32) and the rotation variations (Eq. 33) is 
essential to get a drift on the eccentricity. That is why 
the eccentricity pumping was never observed in previous 
studies that did not take into account the spin evolution. 

The major difference when we consider the full non- 
linearized problem is that the drift (Eq. 37) cannot grow 
indefinitely. Indeed, when the eccentricity reaches high 
values, the drift vanishes (Fig. lb). Moreover, tidal ef- 
fects are also enhanced for high eccentricities and coun- 
terbalance the drift (Eq. 18). As a consequence, the drift 
on the eccentricity is never permanent, although it can 
last for the age of the system (Fig. lc,d). 

In order to observe the pumping effect, the eccentricity 
should not be damped, while the damping timescale of 
the spin of the planet should be of the order of the pe- 



riod of eccentricity oscillation. This is valid for gaseous 
planets roughly within 0.1 < ai < 0.3. About half of 
the planets in this range present eccentricities higher 
than 0.3 (Tab. 1). This is somehow unexpected, since 
the time-scale for damping the eccentricity (r ~ ^/K') 
is shorter than the age of those systems (Tab. 1, Fig. la). 
Thus, unless the initial eccentricity of those planets was 
extremely high, one may suspect of the existence of unde- 
tected companions that help the eccentricity to maintain 
the present values. 

5. CONCLUSION 

In this Letter we have shown that, under some partic- 
ular initial conditions, orbital and spin evolution cannot 
be dissociated. Indeed, counterintuitive behaviors can 
be observed, such as the secular augmentation of the ec- 
centricity. This effect can last over long time-scales and 
may explain the high eccentricities observed for moder- 
ate close-in planets. 

The variations in the flattening of the inner planet due 
to rotation (Eq. 1) is a key element to pump the eccen- 
tricity by means of gx (Eq. 37) . We have considered an 
instantaneous response to the rotation in J2, but a time 
delay until the maximum deformation is reached is to be 
expected. This may increase the phase lag between the 
eccentricity forcing (Eq. 32) and the precession angle os- 
cillations (Eq. 34) . It results that the drift effect on the 
eccentricity can be even more pronounced than the one 
presented here. 
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